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Abstract 

We study the electronic properties (density of states, conductivity and ther- 
mopower) of some nearly-free-electron systems: the liquid alkali metals and two 
liquid alloys, Li-Na and Na-K. The study has been performed within the self- 
consistent second order Renormalized Propagator Perturbation Expansion (RPE) 
for the self-energy. The input ionic pseudopotentials and static correlation func- 
tions are derived from the neutral pseudoatom method and the modified hypernet- 
ted chain theory of liquids, respectively. Reasonable agreement with experiment is 
found for Na, K, Rb and Na-K, whereas for Li and Cs and Li-Na the agreement is 
less satisfactory 
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1 INTRODUCTION 



The introduction of the diffraction model together with the development of the pseu- 
dopotential concept, marked the beginning of extensive calculations of structural and 
electronic properties of liquid metals 0, ^. Based on this concept one can justify the use 
of simple second-order perturbation theory for weak scattering systems like the alkalis. 
The most used expression in this Nearly-Free Electron (NFE) model is the well-known 
Ziman formula for the resistivity. Two key ingredients appear in it, namely the static 
structure factor S{q) and the screened pseudopotential w{q). These two quantities are 
not independent, since the pseudopotential determines the forces between the ions, which 
ultimately determine the structure of the liquid. Therefore S{q) should be obtained from 
w{q) and the result used in the evaluation of the electronic properties, and from the new 
electronic density one should again derive a new pseudopotential. This is what we call a 
self-consistent calculation of the electronic properties of a liquid metal. 

Most of the calculations up to date are, with the exception of ab-initio Molecular 
Dynamics simulations, not self-consistent, i.e. for S{q) one takes either one from an 
analysis of the experimental structure factor data, or one obtained from a model for 
the interatomic interactions, like hard spheres, with parameters fitted to experimental 
data. This model is then used as input either in an approximate scheme to solve the 
Ornstein-Zernike equations like the Percus-Yevick or (Modified) HyperNetted Chain, or 
in a liquid structure computer simulation. Then a pseudopotential is chosen and the 
electronic properties, for example the electronic density of states (EDOS), are calculated. 
There have been some attempts to study the properties in a self-consistent way in a 
more restricted sense, that is, using the same pseudopotential for the structure factor, 
thermodynamic properties and conductivity. However up to the present day such a scheme 
has failed: it leads to disagreement with either the electronic properties, or with the atomic 
properties. 

There are at least two important exceptions in the previous works, which do make 
the calculations within the spirit of self-consistency. One is the work of Jank and 
Hafner , where molecular dynamics are performed with effective pair potentials derived 
from a given pseudopotential, and then some ionic configuration is selected to perform a 
Linearized-Muffin-Tin-Orbitals electronic calculation. 

The other calculations, which are self-consistent, are the modern ab-initio molecular 
dynamics (AIMD) simulations. Using this method, we are only aware, in the case of alkali 
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metals, of EDOS and conductivity calculations for Na and Rb at several temperatures 

In this work, we propose a completely different method, which does not rely on sim- 
ulations, to obtain within a self-consistent scheme, the ionic structure and electronic 
properties of the liquid metal. Moreover, the theory for the calculations of the electronic 
properties which is the second-order Renormalized Propagator Perturbation Expansion 
(RPE) 1^, satisfies the so-called Generalized Optical Theorem (GOT), which provides 
a second level of self-consistency in the theory, namely between the EDOS on the one 
hand and the conductivity on the other hand. This method is much faster and much less 
demanding than AIMD simulations, but, as we will show below, achieves the same level 
of agreement with experiment for weak scattering systems, like the alkali metals, when 
pseudopotentials of similar quality are used. 

A similar RPE-type scheme to calculate the electronic structure has been applied 
before but only in a non self-consistent way (see for a discussion [|], |^, H): the integral 
equation for the self energy was approximated by neglecting the self energy in the right 
hand side of eq. ^]3| below. 

The bare pseudopotentials used in this work were constructed from first principles 
using the Neutral PseudoAtom (NPA) method: one solves a system of one alkali atom 
immersed in an electron jellium in the Local Density Approximation (LDA). Linear Re- 
sponse Theory (LRT) was then applied to obtain the screened pseudopotential w{q) on 
the one hand and the effective pair potential on the other hand; for details see 0, [14| . 
The structure factor S{q) was obtained from the pair potential using the Modified Hy- 
perNetted Chain (MHNC) theory of liquids. This combination has already been applied 
successfully to the study of the static and dynamic structure factors of the liquid alkalis 
as well as the Na-K and Li-Na liquid alloys |T3|, [l^, the systems studied in this paper. 

2 THEORY 

We give here the essentials of the theory, further details can be found in |Q. Atomic units 
with 2m = 1,^ = 1, e^ = l are used throughout. 

We consider a system of ions with their corresponding valence electrons in a volume 
Q. Each valence electron moves in a self-consistent potential due to the ion cores and 
the other valence electrons. Therefore, we can write the single-particle Hamiltonian as: 
H = Ho + V, where Ho is the kinetic energy operator and V(r) = I]jw(|r — Rj|) is an 
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effective one-electron potential which is taken to be the sum of spherically symmetric 
screened local pseudopotentials centered on each ion. The central quantity in this study 
is the one-electron Green function, G{]<.,E). The EDOS per atom, unit energy and unit 
volume is given by 

n(E) = --/^ImG(k,E) (2.1) 

Moreover in terms of the so-called self-energy, S(k, i?), the Green function can also be 
expressed as 

G{k,E) = {E - k"^ -j:{k,E))-^ (2.2) 
Within the second-order RPE |Q, the self-energy is approximated by 

r r/k' 

m,E)=pJ j^y{\k-k'\)S{\k-k'\)G{k',E), (2.3) 

where a constant term pw{q = 0) has been dropped using it as the energy origin, and p 
denotes the atomic number density. This expression for the self energy is easily extended 
to liquid alloys by the substitution of w'^{q)S{q) by J2i,jixiXjY^'^Wi{q)wj{q)Sij{q), where 
Xi and stands for the concentration and Wi{q) for the screened pseudopotential of the 
ith-component, whereas Sij{q) are the Ashcroft-Langreth partial static structure factors. 

Substituting the relation (p.2|) into (|2.3|) , an integral equation is obtained for the self- 
energy, which is solved for each energy, obtaining S(k, £^), G{k,E), and n{E) for that 
energy from equation (|2.1|) . 

The static electrical conductivity is given by the Kubo-Greenwood equation 

a = J dE{-dfFD{E)/dE)a{E) (2.4) 

where fpoiE) is the Fermi-Dirac distribution function and (j{E) is the contribution of 
the electrons of energy E to the conductivity, which is given by 

^^^^ = " ^kAkReG(k, E)) , (2.5) 

where K(k, E) obeys the following integral equation 

Kik, E) ,p /• dk' „ . kk' . , , , 

The vertex function, A(k, k',£'), has still to be specified. A self-consistent calculation of 
the EDOS and the conductivity requires the self-energy and the vertex function to be 
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related by the generalized optical theorem (GOT) 




(2.7) 



We have chosen 



A(k, k', E) = pw^k - k'\)S{\k - k'l) 



(2.8) 



which clearly satisfies the GOT when the second-order RPE expansion is made for the 
self-energy ^]3[ 

Finally, once <j{E) has been obtained, the thermopower, S, can be calculated as a 
function of energy from the Mott relation 



For the calculation of the self energy we used a linear mesh of about 400 points with a 
cut-off of the pseudopotential at ikp- Results of the various electronic properties are well 
converged because a twice as high cut-off does give results which differ only by about 
.1 %. That such a cut-off is sufficient also follows from the observation that the second 
contribution to the conductivity is nearly independent of the energy. 

The rather singular integrand appearing in the expressions for the self energy and the 
conductivity was calculated using a rational expansion up to second order in q in each 
interval of the g-mesh, and were evaluated analytically. Checks using various integration 
schemes showed that this scheme gives a very reliable approximation of these q integrals, 
with an estimated error of .1 % or less. 

The iteration was started from the bottom of the band. Convergence of the self energy 
was obtained within 10 iterations - irrespective of its initialization. No mixing scheme was 
used. The convergence of the integral equation for the conductivity was rather slow near 
the bottom of the band, requiring some 20 to 50 iterations. For the DOS as a function of 
the energy from the bottom of the band to about 2Ef we used about 50 to 100 energy 
points. The solution of the integral equation for the conductivity showed oscillations as 
a function of energy when the g-mesh is taken too coarse. 



S 



7r^k%Td log a{E) 
3|e| Je 



(2.9) 



3 NUMERICAL DETAILS 
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4 RESULTS AND DISCUSSION. 

The preceding formalism has been apphed to hquid alkah metals at thermodynamic con- 
ditions near the melting point as well as to liquid Li-Na and Na-K alloys at the temper- 
atures T = 725 K and T = 373 K respectively. Calculations were also performed for the 
thermodynamic states considered by Silvestrelli et al ^ for liquid Na and Shimojo et al 
for liquid Rb [g. 

The results obtained for the EDOS of the liquid alkalis, and liquid Li-Na and Na-K 
alloys are shown in figure |l[ First, we note that the present results are, in general, very 
similar to those derived from the free electron model and only for Rb and Cs we see some 
structure appearing in the EDOS. The bottom of the band is only slightly shifted below 
the free-electron band; this shift increases on alloying from about 0.003 Ry. for the pure 
elements to 0.02 Ry for the alloys. We checked the validity of these results of the RPE by 
comparing it with results of the Quasi-Crystalline approximation (QCA) and found 
only small changes in the total EDOS with respect to the RPE results. 

The valence bands of the alkalis have been investigated by Indlekofer and Oelhafen 
using various UPS excitation energies. There are two important characteristics observed in 
these valence band spectra: (i) the widths of the valence bands, i. e. the energy difference 
between the Fermi level and the bottom of the valence band, are narrower than the widths 
expected from free-electron behaviour, by some 15 to 25 %, and (ii) the spectra show a 
triangular shape, except for Li where it is parabolic. The same narrowing is observed in 



the crystalline state of alkalis |T^, and has been attributed to electron correlation effects 



not taken into account in the Local Density Approximation |]TT . 

The triangular shape has been shown to be consistent with a free electron EDOS, 
and is explained by the difference in cross sections of the valence electrons with different 
angular momentum for excitations in the UPS regime H. In the RPE these different 
angular moment components are not available. 

Although the present scheme is not able to account for the bandwidth reduction of 
the EDOS, it must be remarked that this is also the case with the AIMD simulations 
reported for liquid Na |^ and Rb In both cases, the valence bands show an almost 
parabolic shape and a width similar to the free-electron model. 



The electrical conductivity is calculated by solving the integral equation |2]^ for -ft'(k, E) 
and using the vertex defined in |2.(j| . The electrical conductivity as a function of energy 
behaves in a similar way for all the systems studied. It is small near the bottom of the 
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band, and then it increases to a broad maximum around the Fermi level, and then it de- 
creases. The Fermi level is in all cases located before the maximum in (j{E), so according 
to relation |2.9| leading always to a negative thermopower (see table |l|). The sign of the 
thermopower agrees with the experimental values except for Li and Cs. 

In table |I| we present the resistivity we calculated in the present RPE scheme together 
with those predicted by Ziman's formula using the same structure factor and pseudopoten- 
tial. The values of the latter approach are in most cases somewhat smaller than the RPE. 
The comparison with experimental data reveals discrepancies for Li and Cs, however, the 
discrepancies are smaller in the case of alloys. 

In the case of Na at different temperatures we find that, although the structure factor 
is well reproduced, the values obtained for the conductivity are always smaller than the 
experimental ones, the discrepancy with experiment becoming smaller as the temperature 
increases. The same behaviour for the resistivity is obtained by Silvestrelli et al 0] in 
their AIMD studies for Na; moreover, the discrepancy with experiment is in their AIMD 
calculation larger for lower temperatures than in the RPE. They used the Topp-Hopfield 
pseudopotential for Na, which reproduced correctly the structure factor in all thermo- 
dynamic states. Moreover they performed computations with norm-conserving non-local 
pseudopotentials, even including core corrections, but no improvement was obtained for 
the conductivity values, compared with experiment. 

For Rb, the AIMD of Shimojo et a/ [Q, using a norm-conserving non-local TrouUier- 
Martins pseudopotential, obtained conductivities which are much smaller than the ex- 
perimental ones, whereas our calculations agree very well with experiment. It is also 
interesting to note that the structure factor was correctly reproduced both in the AIMD 
and in our calculations. 

We also performed calculations for the atomic and electronic properties using the 
Fiolhais-Perdew pseudopotential but did not obtain improved agreement with ex- 
periment . 

5 CONCLUSIONS. 

The above comparisons show that when suitable pseudopotentials are used, the com- 
bination of LRT, MHNC and RPE produces self-consistent results for the atomic and 
electronic properties of liquid alkali metals, which agree with modern state of the art 
AIMD methods, at least for alkalis for which results using both these methods are avail- 
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able. The agreement of our atomic structure factors with experiment is excellent. The 
agreement between theoretical and experimental electronic properties is less satisfactory. 
For Li and Cs, we find large discrepancies with experimental data for the resistivity and 
for thermopower of Li and Cs for we do not even find the right sign. No thermopower 
data are available from these AIMD simulations. 

In our calculations we have found that for the pure alkali metals, the deviation of 
the EDOS from the free electron parabola is rather small, and it increases somewhat on 
alloying. The resistivity and thermopower agree reasonably well with experiment for Na, 
K, Rb and Na-K but is less satisfactory for Li and Cs, and as a consequence for the liquid 
Li-Na alloys in the Li-rich composition regime. It is known that the pseudopotential 
of Li and Cs has rather large non-local components, for Li because p states are missing 
in the core, and for Cs because d states are present near the Fermi level. Such non- 
local components of the pseudopotential do not seem to contribute much to the atomic 
structure, but could be of importance for the electronic properties. However, the use 
of a non-local pseudopotential by Shimojo |^ does not improve the AIMD result for 
the conductivity for Rb when compared with experiment. Note that for Rb we find 
good agreement with experiment. Such non-local components could alter the position 
of the maximum in a as a function of energy. Introducing relative stronger electron-ion 
scattering below the Fermi level, which would give a minimum in cr{E), thus also changing 
the sign of the thermopower. 

The larger bandwidth calculated for the EDOS compared with experiment has been 
attributed to electron correlations not taken into account in the LDA approximation. 
In our scheme we take these electron-electron interactions into account implicitly in the 
screened pseudopotential. To take these electron-electron interactions into account in a 
more explicit way would require the introduction of a self-energy in the r.h.s. of equation 



Ol , and using the bare electron-ion potential: in terms of diagrams it means that instead 
of summing part of the electron-electron interactions in a screened electron-ion potential, 
one collects them in a medium propagator, which describes the motion of electrons includ- 
ing the electron-electron interactions between two scattering events on bare electron-ion 
potentials. This electron self-energy can be approximated for these NFE systems in the 
Random Phase Approximation (see [0). The consequence of such a scheme for the con- 
ductivity is not an increase of the effective electron-ion scattering, but could result in a 
slight change of the EDOS at the Fermi level. Most of the changes in the EDOS should 
occur below the Fermi level. The EDOS at the Fermi level of Na is nearly given by the 
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NFE model. 

The discrepancy of the electronic transport properties with experimental results we 
attribute to the construction of the pseudopotential, and not to basic problems in our 
theory which should be valid for these weak scattering systems, as also follows from the 
comparison of our results for the resistivity and those obtained from the Ziman equation: 
we find only small corrections to the Ziman results. 

The success of both the present approach and the AIMD in reproducing the exper- 
imental electronic properties appears to depend very much on the subtle details of the 
pseudopotential. The atomic structure and electronic structure and transport properties 
are determined by different parts of the pseudopotential. Electrons penetrate more into 
the core region and are also more sensitive to angular non-local components. So changes 
in these contributions not necessarily will affect the atomic structure. 

The conclusion is that there exists accurate theoretical tools that can in principle 
give reliable self-consistent results for both the atomic and electronic properties of liq- 
uid (alkali) metals provided adequate pseudopotentials are used. However, considering 
the results of the present calculations, we conclude that most of the pseudopotentials 
proposed up to now for alkali metals, although good enough to describe the static and 
dynamic atomic structure, are not accurate enough to give good agreement also with the 
experimental electronic properties. Even for Na, which has always been considered as the 
simplest metal to describe. 

In contrast to the pseudopotentials used in AIMD simulations our pseudopotentials 
arc adjusted to the electron density in the liquid metal. We attribute the failure of our 
pseudopotentials to describe the EDOS to the linear screening approximation. Further 
checks are currently being undertaken. Preliminary results for a pseudopotential based 
on the electronic density show better agreement with experiment for the electronic prop- 
erties, while preserving the good agreement of the atomic properties with experiment. 
We attribute the discrepancy between experiment and theory in case of the electronic 
transport properties to so-called nonlocal contributions to the pseudopotential, which in 
case of Li and Cs can be rather large. 

We have shown that fully self-consistent ab-initio calculations, with as only input the 
chosen pseudopotential, using integral equation approximations for the static and dynamic 
atomic structure as well as for the electronic properties give results of the same reliability 
as AIMD simulations, under the condition that the electron-ion pseudopotential has been 
constructed well. The advantage of the present method is that it is much faster than the 
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AIMD simulations and only requires minimal computational resources, and so the study of 
possible improvements in the construction of the pseudopotential taking into account non- 
local components and implementation of a scheme bare electron-ion potentials together 
with the electron-electron self energy can in principle easily be implemented and studied. 
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Figure 1: The EDOS as a function of energy for the liquid alkalis near the melting point, 
and for the hquid alloys Li-Na at T=725 K and Na-K at T=373 K. 
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Table 1: Calculated resistivity (nflcm) and thermopower ( fiV/K) for the liquid alkali 
metals near melting, the liquid alloys Li-Na at T—72b K and Na-K at T=373 K. 





p(Kubo) 


p(Ziman) 


p(exp) 


S(Kubo) 


S(exp) 


Li 


7.3 


7.0 


24 


-3.9 


+21.7 


Na 


16.3 


15.8 


9.6 


-8.0 


-7.9 


K 


19.6 


18.5 


13.0 


-10.6 


-14.0 


Rb 


22.0 


20.8 


22.0 


-11.2 


-6.3 


Cs 


14.9 


13.8 


36.0 


-11.7 


+6.4 


Na 


36.2 


36.4 


24.6 


-16.0 


-13.2 


Lio.4Nao.6 


27.8 


27.5 


30.2 


-13.9 


+2.2 


Lio.sNao.s 


25.9 


25.3 


31.5 


-13.4 


+5.1 


Lio.6Nao.4 


23.2 


22.8 


32.2 


-12.6 


+8.5 


Lio.8Nao.2 


18.2 


17.8 


33.1 


-11.1 


+ 16.2 


Li 


12.3 


12.0 


33.6 


-9.0 


+25.2 


K 


21.6 


20.3 


14.9 


-11.5 




Nao.sKo.T 


48.4 


48.2 


39.5 


-8.5 




Nao.5Ko.5 


48.4 


49.0 


43.0 


-8.6 




Nao.8Ko.2 


32.6 


32.8 


26.0 


-8.3 




Na 


16.3 


15.8 


9.6 


-8.0 
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